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ABSTRACT 

We use a statistical thermodynamic approach to determine the composition of clathrate hydrates which may form from 
a multiple compound gas whose composition is similar to that of Titan's atmosphere. Assuming that noble gases are 
initially present in this gas phase, we calculate the ratios of xenon, krypton and argon to species trapped in clathrate 
hydrates. We find that these ratios calculated for xenon and krypton are several orders of magnitude higher than in 
the coexisting gas at temperature and pressure conditions close to those of Titan's present atmosphere at ground level. 
Furthermore we show that, by contrast, argon is poorly trapped in these ices. This trapping mechanism implies that 
the gas-phase is progressively depleted in xenon and krypton when the coexisting clathrate hydrates form whereas the 
initial abundance of argon remains almost constant. Our results are thus compatible with the deficiency of Titan's 
atmosphere in xenon and krypton measured by the Huygens probe during its descent on January 14, 2005. However, 
in order to interpret the subsolar abundance of primordial Ar also revealed by Huygens, other processes that occurred 
either during the formation of Titan or during its evolution must be also invoked. 
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' 1. Introduction 

. An unexpected feature of the atmosphere of Titan is that 

■ no primordial noble gases other than argon were detected 
by the Gas Chromatograph Mass Spectrometer (GCMS) 

. aboard the Huygens probe during its descent on January 

■ 14, 2005. The detected argon includes primordial ^^Ar (the 
' main isotope) and the radiogenic isotope '^''Ar, which is a 

decay product of (Niemann et al. 2005). The other 
primordial noble gases ^^Ar, Kr and Xe were not detected 
, by the GCMS instrument, yielding upper limits of 10~^ 

■ for their mole fractions. In any case, •^^Ar/^''N is about six 
' orders of magnitude lower than the solar value, indicating 

that the amount of ^^Ar is surprisingly poor within Titan's 
atmosphere (Niemann et al. 2005). 
, Several scenarios have been proposed in the literature 
to interpret the deficiency of Titan's atmosphere in primor- 
dial noble gases. The first category of scenarios postulates 
that this deficiency results from processes that occurred ei- 
ther during Titan's accretion or during the formation of its 
planetesimals. In this way, Alibert & Mousis (2007) pro- 
posed that the building blocks of Titan initially formed 
at low temperature in the solar nebula and were subse- 
quently partially vaporized in Saturn's subnebula during 
their migration and accretion, thus leading to the loss of 
noble gases and of carbon monoxidfQ. Alternatively, Owen 
(2006) proposed that Titan's planetesimals were directly 
produced at temperatures high enough in Saturn's subneb- 
ula to impede the trapping of the noble gases during their 
formation and growth. The second category of scenarios 
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^ CO is several orders of magnitude less abundant than CH4 
in Titan's atmosphere. 



postulates that Titan initially incorporated primordial no- 
ble gases during its formation but that subsequent processes 
in its atmosphere or on its surface prevented them from be- 
ing detected during the descent of the Huygens probe. For 
example, Jacovi et al. (2005) suggested that the aerosols 
observed in Titan's atmosphere may have cleared its con- 
tent in noble gases, assuming they were produced from the 
aggregation of polymers. 

Another interpretation is the trapping of atmospheric 
compounds in clathrate hydrates (Osegovic & Max 2005 
- hereafter OM05). Indeed, in such ice structures, water 
molecules form cages which are stabilized by the trapping 
of volatiles. Thus, OM05 calculated clathrate hydrate com- 
positions on the surface of Titan by using the program 
CSMHYD developed by Sloan (1998) and showed that such 
crystalline ice structures may act as a sink for Xe. However, 
from these results and because no other noble gases are con- 
sidered in this program, they assumed similar trapping effi- 
ciencies for Ar and Kr. In addition, because the CSMHYD 
code is not suitable below 160 K for gas mixtures rele- 
vant to the atmospheric composition of Titan, OM05 also 
extrapolated their results to the surface temperature of 
Titan (about 94 K). Moreover, the stability curves of multi- 
ple guest clathrate hydrates calculated from the CSMHYD 
code strongly depart from the experimental data at low 
temperature. 

Here, we reinvestigate OM05's assumptions and calcu- 
late more accurately the trapping of noble gases in clathrate 
hydrates using a statistical thermodynamic model based 
on experimental data and on the original work of van der 
Waals & Platteeuw (1959). In particular, the relative abun- 
dances of Xe, Kr, and Ar are explicitly calculated in condi- 
tions that are valid for the present temperature and pres- 
sure at the surface of Titan. Our results show that, when 
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clathration occurs at low temperature and pressure, the 
trapping of Xe and Kr is very efficient, contrary to that of 
Ar. This conclusion is partly consistent with the Huygens 
probe measurements, and supports indirectly the presence 
of noble gas-rich clathrate hydrates on the surface of Titan. 

2. Clathrate hydrates composition 

We follow the method described by Lunine & Stevenson 
(1985) to calculate the abundance of guests incorporated 
in clathrate hydrates from a coexisting gas of specified 
composition. This method uses classical statistical mechan- 
ics to relate the macroscopic thermodynamic properties of 
clathrate hydrates to the molecular structure and interac- 
tion energies. It is based on the original ideas of van der 
Waals & Platteeuw (1959) for clathrate formation, which 
assume that trapping of guest molecules into cages corre- 
sponds to the three-dimensional generalization of ideal lo- 
calized adsorption. In this model, the occupancy fraction of 
a guest molecule K for a given type t of cage {t = small or 
large) in a given type of clathrate hydrate structure (I or 
II) can be written as 



yK,t 



Ck.iPk 



(1) 



where the sum at the denominator includes all the species 
which are present in the initial gas phase. CK,t is the 
Langmuir constant of species K in the cage of type t, and 
Pk is the partial pressure of species K. This partial pres- 
sure is given by Pk = xk x P, with xk the molar fraction of 
species K in the initial gas phase, and P the total pressure. 

The Langmuir constants are determined by integrating 
the molecular potential within the cavity as 
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r dr. 
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where Rc represents the radius of the cavity assumed to 
be spherical, and w{r) is the spherically averaged Kihara 
potential representing the interactions between the guest 
molecules and the II2O molecules forming the surrounding 
cage. Following McKoy & Sinanoglu (1963), this potential 
w{r) can be written for a spherical guest molecule, as 
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The Kihara parameters a, <j and e for the molecule-water 
interactions of the species considered here are given in Table 
[TJ In Eq. ([3]), 2; is the coordination number of the cell and 
r the distance of the guest molecule from the cavity center. 
These parameters are given in Table [5] and depend on the 
structure of the clathrate and on the type of the cage (small 
or large). Note that the Kihara parameters are based on fit 
to experimental data. However, a systematic study of the 



Table 1. Parameters of Kihara potential: a is the radius of 
the impenetrable core, e is the depth of the potential well, 
and a is the Lennard-Jones diameter. 



Molecule 


-(A) 


a(A) 




Ref. 


CH4 


3.1514 


0.3834 


154.88 


(a) 


C2H6 


3.2422 


0.5651 


189.08 


(a) 


N2 


3.0224 


0.3526 


127.67 


(a) 


Ar 


2.829 


0.226 


155.30 


(b) 


Kr 


3.094 


0.307 


212.70 


(b) 


Xe 


3.3215 


0.2357 


192.95 


(a) 



(a) Jager (2001) 

(b) Diaz Pena et al. (1982) 

Table 2. Parameters for the cavities, h represents the num- 
ber of small ihs) or large (6^) cages per unit cell for a given 
structure of clathrate, Rc is the radius of the cavity (param- 
eters taken from Sparks et al. 1999), and z the coordination 
number in a cavity. 



Clathrate structure 




I 




II 


Cavity type 


small 


large 


small 


large 


b 


2 


6 


16 


8 


Rc{k) 


3.905 


4.326 


3.902 


4,682 


z 


20 


24 


20 


28 



influence of their values on the present results would be 
required to assess the relevance of the present conclusions 
(Papadimitriou et al. 2006). 

Finally, the relative abundance fpc of a guest molecule 
K in a, clathrate hydrate can be calculated with respect to 
the whole set of species considered in the system as 



bsUK.s + biyK,i 



(5) 



where hs and are the number of small and large cages per 
unit cell respectively, for the considered clathrate hydrate 
structure (I or II). The sums in the denominator include 
all species present in the system. Note that, even if it has 
been shown that the size of the cages can be reduced by ~ 
10-15% at low temperatures (Shpakov et al. 1998; Besludov 
et al. 2002), we assume here that this size remains fixed. 
In the same way, we neglect possible cage distortion due to 
large guests. Moreover, the relative abundances calculated 
this way do not discriminate between the various isotopes 
of the same species. 

All calculations have been performed at the dissocia- 
tion pressure P = P^j^f of multiple guest clathrate hy- 
drate, i.e. temperature and pressure conditions at which 
the clathrate hydrate is formed. This dissociation pressure 
can be deduced from the dissociation pressure Pf^^^ of a 
pure clathrate hydrate of species as (Hand et al. 2006) 



(4) 



diss 



\ / ^ pdiss / 



(6) 



The dissociation pressures P^^^^ derived from laboratory 
measurements follow an Arrhenius law (Miller 1961) as 



logiP^'n 



(7) 



where P'^^^^ is expressed in Pa and T is the temperature in 
K. The constants A and B fit to experimental data (Lunine 
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& Stevenson 1985; Sloan 1998) in the present study are 
given in Table [3] Note that, in the present approach, the 
dissociation pressures and the Langmuir constants are in- 
dependently calculated. 

Table 3. Parameters of the dissociation curves for various 
single guest clathrate hydrates. A is dimensionless and B is 
in K. 



Molecule 


A 


B 


CH4 


9.89 


-951.27 


C2H6 


10.65 


-1357.42 


N2 


9.86 


-728.58 


Ar 


9.34 


-648.79 


Kr 


9.03 


-793.72 


Xe 


9.55 


-1208.03 



3. Results and discussion 

We determine the relative abundances /xo, /Kr and /Ar 
with respect to all guests incorporated in clathrate hydrate 
formed from a gas whose composition is similar to that of 
Titan's atmosphere. Similarly to OM05, we assume that 
N2, CH4, C2H6 and the considered noble gas (Xe, Kr or 
Ar) are the main atmospheric compounds of the initial 
gas phase (see Table H]). Note that when the calculations 
are performed for one noble gas (Xe, Kr or Ar), the two 
others are excluded from the initial gas composition. The 
gas phase abundance of N2 is taken similar to OM05 and 
that of CH4 has been updated from the Huygens measure- 
ments (Niemann et al. 2005). The atmospheric abundance 
of ethane has been fixed to 0.1% in our calculations. For 
each considered noble gas, three initial abundances have 
been selected in order to investigate the influence of its 
initial concentration in the gas phase, on the clathrate hy- 
drate's trapping efficiency. The abundance of CH4 is then 
varied in consequence to preserve the global volatile bud- 
get. Moreover, we have calculated the relative abundances 
of noble gases for both clathrate hydrates structures I and 
IL Because the conclusions of our calculations are similar 
for both structures, we focus here on the results obtained 
for structure I only. Note that this structure corresponds to 
that derived by OM05 from the CSMHYD program. 

Table 4. Initial gas phase abundances considered in the 
atmosphere of Titan. For each noble gas (Xe, Kr and Ar), 
three different compositions of the initial gas phase are in- 
vestigated. 



Molecule Molar fractions (%) 

Noble gas OJ 0055 OOT 

CH4 4.8 4.845 4.89 

N2 95 95 95 

C2H6 0.1 0.1 0.1 



Figures 1 to 3 represent the evolution of /xc, /Kr and 
/Ar respectively, in a multiple guest clathrate hydrate as 
a function of T along the dissociation curve P^f^ = f{T), 
and for three different noble gas abundances in the initial 
gas phase (see Tabled). Figure [T] shows that, irrespective 
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Fig. 1. fxe hi a multiple guest clathrate hydrate, as a func- 
tion of temperature and of the initial gas phase composition 
(see Table HI). 
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Fig. 2. Same as FigUbut for Kr. 



of the initial gas phase abundance of Xe, its relative abun- 
dance progressively increases in clathrate hydrate when T, 
and hence P^'^f , decrease. This behavior of the Xe clathra- 
tion cfhciency agrees with the calculations of OM05 (e.g. 
their Fig. 5). Similar behavior is seen with Kr, as shown 
in Figure [21 By contrast, the relative abundance of Ar in 
clathrate hydrate decreases when T and P^^f^ decrease, as 
shown in Fig. [3l This figure also shows that Jat gradually 
converges towards zero at temperatures below ~ 100 K, ir- 
respective of its initial concentration in the gas phase. This 
huge difference between the trapping of Ar, Xe and Kr in 
clathrate hydrate is mostly due to the higher values of the 
Kihara parameters of the two latter molecules (see Table 
[1]) which are responsible for a stronger interaction energy 
with the surrounding water molecules in clathrate hydrate 
cages. 

In addition to the relative abundances discussed 
above, we have also calculated the abundance ratio for the 
three considered noble gases Xe, Kr and Ar (Table [5|). This 
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Fig. 3. Same as Fig[T]but for Ar. 

Table 5. Abundance ratios of noble gas in clathrate hy- 
drate to noble gas in the initial gas phase for Xe, Kr and 
Ar. These ratios are calculated at r=171 K and P^f^ = 
1.5 bar. 



Initial molar fraction 


abundance 


in gas 


ratio 


Xe gas 




0.001 


168 


0.00055 


184 


0.0001 


204 


Kr gas 




0.001 


195 


0.00055 


285 


0.0001 


555 


Ar gas 




0.001 


1.577x10"^ 


0.00055 


1.572x10"^ 


0.0001 


1.566x10"^ 



ratio is defined as the ratio between the relative abundance 
/k of a given noble gas in the multiple guest clathrate hy- 
drate and its initial gas phase abundance xk (see Table H]). 
We considered one particular point located on the dissoci- 
ation curve of the multiple guest clathrate hydrate under 
study in our calculations. This point corresponds to the 
present atmospheric pressure at the ground level of Titan 
(P^ix = 1-5 bar) and to a temperature value of 171 K, as 
indicated by the clathrate hydrate dissociation curve. Table 
[5] shows that, depending on the initial gas phase conditions 
and for this particular point, the relative abundances of Xe 
and Kr trapped in clathrate hydrates are between 168 and 
555 times those in the initial gas phase. Considering the 
current amount of atmospheric CH4 calculated by Lunine 
& Tittemore (1993) (~ 3 x lO^" g), we have calculated that 
a 0.1% gas phase abundance of a given noble gas trans- 
lates into the presence of - 3.9 X lO^"^ moles in Titan's 
atmosphere. Assuming full efficiency for the sink mecha- 
nism described here, the amount of water needed for the 
trapping of Xe or Kr in clathrate hydrates at the surface 
of the satellite is ~ 4 x 10^^ g. This roughly corresponds to 



a single guest (Xe or Kr) clathrate hydrate layer ~ 50 cm 
thick at the surface of Titan. 

Such an efficient trapping may be high enough to signif- 
icantly decrease the atmospheric concentrations of Xe and 
Kr, provided that clathrate hydrates are abundant enough 
on the surface of Titan. On the contrary, with an abundance 
ratio much lower than 1, Ar is poorly trapped in clathrate 
hydrate and the Ar atmospheric abundance consequently 
remains almost constant. 

4. Conclusions 

The trapping mechanism we investigated can explain the 
deficiency in primordial Xe and Kr observed by Huygens in 
Titan's atmosphere. Indeed, our results show that these 
noble gases could have been progressively absorbed in 
clathrate hydrates located at the surface of Titan during 
its thermal history. By contrast, our results show that the 
trapping of Ar is poor in clathrate hydrates. As a conse- 
quence, such a mechanism alone cannot interpret the sub- 
solar abundance of primordial Ar which is also observed 
in Titan's atmosphere. The different scenarios invoked in 
the Introduction may thus have worked together with the 
trapping in clathrate hydrates to explain the deficiency in 
primordial Xe, Kr and Ar. Note that boulders observed 
in images from Huygens have been invoked as evidence 
for the existence of clathrate hydrates on the surface of 
Titan (OM05). However, further in situ measurements are 
required to investigate whether noble gases-rich clathrate 
hydrates really exist on the surface of Titan. 
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